Flow Equations for the Ionic Hubbard Model 



oo 

o 

O 



O 
CO 



c3 



i 

c 
o 
o 



> 

CO 

oo 

o 

oo 
o 



% 



Mohsen Hafez, 1 S. A. Jafari*, 2 ' 3 and M. R. Abolhassani 1 

'Department of Physics, Tarbiat Modares University, Tehran, Iran 

'Department of Physics, Isfahan University of Technology, Isfahan 84154-83111, Iran 

3 The Abdus Salam ICTP, 34100 Trieste, Italy 

Taking the site-diagonal terms of the one-dimensional ionic Hubbard model (IHM) as Ho, we employ Con- 
tinuous Unitary Transformations (CUT) to obtain a "classical" effective Hamiltonian in which hopping term 
has been integrated out. For this Hamiltonian spin gap and charge gap are calculated at half-filling and subject 
to periodic boundary conditions. Our calculations indicate two transition points. In fixed A, as U increases 
from zero, there is a region in which both spin gap and charge gap are positive and identical; characteristic of 
band insulators. Upon further increasing U, first transition occurs at U = U C1 , where spin and charge gaps 
both vanish and remain zero up to U — U C2 ■ A gap-less state in charge and spin sectors characterizes a metal. 
For U > U C2 spin gap remains zero and charge gap becomes positive. This third region corresponds to a Mott 
insulator in which charge excitations are gaped, while spin excitations remain gap-less. 

PACS numbers: 71.10.Fd,78.66.Nk,77.22.Ej 



INTRODUCTION 

Ionic Hubbard model (IHM) has been used to study neutral- 
ionic transition in organic compounds JjJ] and understanding 
role of strong correlations in ferroelectricity of metal oxides 
such as BaTiC>3, KNbC>3, KTaC<3 |2|], as well as some quasi 
one-dimensional materials such as (TaSe^I, K0.3M0O3 y|]- 
IHM Hamiltonian includes an staggered one-body external 
potential in addition to the Hubbard Hamiltonian. The Hamil- 
tonian for this model is as follows: 

h = -t^2(c\ a c i+1<7 + h.c.) + uj24Ai ci L c ^ 

ia i 

A 



^(-lYclcir, 



(1) 



where c. i<T (c ia ) is the usual annihilation (creation) operator 
at site i with spin a, t the nearest neighbor hopping ampli- 
tude, U the on-site coulomb interaction parameter, and A a 
one-body ionic potential. This model has a very interesting 
phase diagram at zero temperature at half-filling. When A <C 
U this Hamiltonian, like the usual Hubbard Hamiltonian, is 
transformed to the Heisenberg spin Hamiltonian [l|, |5[] that 
describes a Mott insulator. In the opposite limit A ^> U, 
one can ignore U, and a Bogoliubov transformation gives a 
simple band insulator with gap A. The limit U ~^> A is a 
many-body insulator, while the limit U <C A is a one-body 
insulator. Hence this model is a basis to study the issues like 
matrix-element effects in optical spectra of band insulators, 
versus many -body insulators jgj. 

Apart from offering two entirely different insulating states, 
yet there remains interesting question of the intermediate 
phase: What is the nature of the ground state for U ~ A? 
Researchers have used various methods such as exact diago- 
nalization (ED) yl|4J], density matrix renormalization group 
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(TJMRG) A S H, quantum Monte Carlo (OMQ fil [ll El 
ll], dynamical mean field theory (DMFT) IU3J[l4[l5ll. etc. to 
study the properties of this model. 

The nature of intermediate phase still remains controver- 
sial: Bosonization study of Fabrizio and coworkers 11611 indi- 
cates an spontaneously dimerized phase, which is also sup- 
ported by spin-particle transformation study of Batista and 
coworker s II 1711 . There are also other works supporting this 
scenario ITi [a [ill [l5[]. On the other hand, there are studies 
showing that the intermediate phase is metallic II lL |4L 1 1 OL 1 1 3 1 - 
Brune finds a metallic transition point in between the band 
and Mott insulating states, with simultaneous bond order |9|]. 
Craco et. al. in addition to metallic region, report on a coexis- 
tence phase between the band and Mott insulating states 11411 . 

In this work, we employ the method of continuous uni- 
tary transformations or flow equations to study the half-filled 
one-dimensional IHM in zero temperature subject to periodic 
boundary conditions. CUT method is used to obtain a "clas- 
sical" effective Hamiltonian for IHM. This effective Hamil- 
tonian which is free of fermionic minus sign problem, pro- 
vides a simple picture of the nature of excitations when in- 
terpreted in terms of the transformed ground state. Then for 
L = 20 sites and with A = 20 (energies are in units of 
t) spin and charge gap are numerically calculated for differ- 
ent values of U. These calculations indicate two transition 
points. For U < U Cl , spin and charge gap are both positive, 
for U C1 < U < U C2 are both zero and for U > U C2 spin gap 
remains zero while charge gap becomes positive. Hence this 
the transformed effective Hamiltonian gives band insulator for 
small U region (U < U Cl ), Mott insulator for large U regime 
(U > U C2 ), and metal for the region in between. 



ATOMIC LIMIT 

In the limits U <C A and Act/, the model becomes 
easy to understand and gives band and Mott insulators, respec- 
tively. Yet, another interesting limit in which the phase dia- 
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FIG. 1 : (Color online) Lowest energy configurations for (a) Perfect 
band insulator, (b) Perfect Mott insulator. 



gram becomes particularly simple is the atomic limit (t = 0), 



A 



E[n la ] = -J2(-V l n l „ + Uj2 



riiinii, 



(2) 



where the Hamiltonian can be written solely in terms of classi- 
cal (commuting) variables m a . In this sense, the atomic limit 
corresponds to the classic limit of a quantum problem, and it 
becomes much easier to analyze. Any configuration \{rii a }) 
is an eigen-state of this Hamiltonian. The hopping term which 
causes quantum fluctuations of various |{nicr}) configurations 
by admixing them is absent here. Hence the problem reduces 
to finding out the configuration with lowest energy. To this 
end, we consider a unit cell in real space of the ID lattice 
shown in Fig. Q] Since this Hamiltonian contains purely lo- 
cal degrees of freedom, we consider a zero momentum mode 
(k = 0) which corresponds to repeating this unit cell all over 
the ID chain. In part (a), a state with energy E a = U — A 
per unit cell is shown, while the one in part (b) has energy 
Eb = per unit cell. For E a < Eb, i.e. U — A < 0, state (a) 
indicating a band insulator is stabilized, while for U > A, the 
Mott insulating state with energy Eb is stabilized. At U = A 
these two states are degenerate, and hence in k = mode, j 
spin can freely move in the unit cell of Fig. Q] Therefore the 
picture in the atomic limit is as follows: Eq. (O describes a 
Mott insulator for A < U and a band insulator for A > U. 
The transition line A = U is a metal phase shrinked to a line. 
The question we would like to address in this paper is, what 
happens if we turn on the quantum fluctuations around the 
Hamiltonian (J2]i by introducing a finite hopping t? It turns out 
that CUT is a very suitable tool to investigate this strong cou- 
pling problem. We advertise the result in advance: The inclu- 
sion of the hopping t will renormalize U — ► U(t, U, A), A — > 
A(i, U, A), t — > t — 0, as well as generating new couplings. 



CONTINUOUS UNITARY TRANSFORMATIONS 



18] was intro- 
19D, and Weg- 



Flow equations approach or, CUT method 
duced independently by Glazek and Wilson | 
ner M20I1 . Since then, this method has been applied extensively 
to study various models in condensed matter and high energy 
physics. Some examples include, the Anderson model 12 ill . 



Lipkin model 022I1 . mapping of the electron-phonon interac- 
tion to an effective electron-electron interaction [230, Sine- 
Gordon model II24I1 . RKKY interaction ||25|| . as well as non- 
equilibrium dynamics of strongly correlated systems II26I1 and 
an attempt for a unified description of Fermi and Luttinger 
liquids in all energy scales 02711 . 

In this method the Hamiltonian is transformed by a unitary 
operator U{£) where £ is a parameter between zero to infinity. 
Flow equation for H{£) = W{£)HU{£) becomes, 



d i H(£) = [r ) (£),H(£)}, 



(3) 



where r)(£) = d(U^{£).U{£) is an anti-Hermitian operator 
called generator. Here U(£ = 0) = 1, as at £ = transformed 
Hamiltonian is equal to H . Wegner suggested the following 
generator [20] : 



r,(£) = [H d (£),H r (£)}, 



(4) 



where H d {£) (H r (£)) is diagonal (off-diagonal) part of the 
Hamiltonian. With this choice for the generator, transformed 
Hamiltonian flows towards a diagonal or block-diagonal form 
in the £ — ► oo limit. Mielke introduced another generator that 
is specially useful for band matrices IJ28I1 and preserves the 
band nature of the Hamiltonian matrix. In this section, we use 
Wegner generator to obtain an effective Hamiltonian for IHM. 
It is obvious from Eq. (O that the flow equations without 
any approximation is very hard to solve for Hamiltonians that 
include interaction. In applying the flow equations to IHM we 
use an approximation similar to the one used by Kehrein and 
Mielke in the context of Anderson model J2lll . We approxi- 
mate H{£) as, 



H{£) = ~t(£)J2(clc l+ ia + h. 



Mi) 



e^^u 



+My c t c t 



iu/CiaiCicr 



-V(£) E c \A+\at C i+l<r,Cio 



(5) 



with initial conditions t(0) = 1, A(0) = A, C/(0) = U 
and V(Q) = 0. Physical motivation behind considering n.n. 
Coulomb interaction V is as follows ||8|] : In the simple phase 
diagram shown in Fig.Q]for a fixed A, as one increases U, the 
ionicity ns — ua which is the population difference between 
A, B lattice sites, suddenly drops from 1 to at U C2 = A. If 
an attractive n.n. term V can be generated, it opens up an in- 
termediate phase in which the < ub — nA < 1; i.e. the unit 
cell gently looses its polarization with respect to the half-filled 
background and turns into a Mott insulator. 

In Eq. (O, the first term is taken to be off diagonal and other 
terms are diagonal. This choice resembles the strong coupling 
perturbation schemes starting from the atomic limit. There 
one usually has to organize the hopping processes into those 
terms which change the double occupancy, and those which 
do not [5]. But here we show that the inclusion of smallest 



amount of A makes it possible to renormalize the whole hop- 
ping processes to zero (Eq.fPfli. 

Using equation (@), With this convention, the Wegner gen- 
erator is obtained as: 

r,{£) = t(£)A(£)J2(-lY(4 + i^ - h.c.) 

icr 

-t{£)U{£) ^(cJ <7 cJ .,c io . / Ci + i (r - c\_ la c\ al c l(Jl c l(J - h.c.) 

icrcrt 

-t{£)V(£) ^{c\ a c\ +lrTl C l+lcrl C l+la + 4 4T ct +1<r/ Ci+2<rlCuT 

iaa' 

- C ia C ial C i+^' C i<? _ c i _ lo .c i+lo . / Ci+i CT /Cj . - h.c.) (6) 

To obtain the flow equations, we need to calculate the com- 
mutator of r/(£) and H(£). To do so, we split ?/(£) into three 
parts 



v(£) = mW + m(e) + m(£), 



(7) 



where rji(£), rj2(£), and 773 (£) denote the first, second, and 
third terms in equation (O respectively. Similarly, H (£) in 
equation © has four terms 

H{£) = H X {1) + H 2 {£) + H 3 (£) + H^t), (8) 

where Hi(£) . . . H±{£) correspond to four terms in Eq. ©. 
Therefore: 

W),H(£)\ = [ m (£),H 1 (£) + H 2 (£)} + [ m (£),H 1 (£)} 
+ [773 (£) , Hi {£)] + irrelevant terms, (9) 

where "irrelevant terms" are newly generated couplings which 
are not similar to Eq. (0). The commutators are calculated as 
follows: 

[rj(£),H 1 (i) + H 2 (£)] = 

2t 2 (£)A(£) ]T(-l)*( c t CTC , ia + c l Cl+2a + h. c .) 



and 



+t(£)A 2 (£)Y,(c\ +la c l „ + h.c.), 



teW.^Wh^Mf/WE^E^ 



(10) 



Substituting in Eq. ( ITTb and dropping irrelevant terms gives, 
[r l2 (£),H 1 (£)]=2t 2 (£)U(£) (12) 

x ^(clAa'Cio'Ci-y ~ 4+ia c \a> c i<y' c i+i<y + h - c -) 

iaa' 

The last commutator up to irrelevant terms becomes: 

[ m (£),H 1 (£)] = 2t 2 (£)V(£) (13) 

x *}2{ 2c \A+i<y> c i+i<y' c i<y - c\ a c\ a ,c la/Cia + h.c.) 

iaa' 

Substitution of Eqs. ( fTUt . (fT2l) and ( fT3b in Eq. ©, gives the 
following set of flow equations 

d e t(£) = -t(£)A 2 (£) (14) 

d e A(£)=8t 2 (£)A(£) (15) 

d e U(£) = 8t 2 (£) (U(£) ~V{£)) (16) 

(17) 



d e V(£) = 4t 2 (£){2V{£) -U{t)) 



We obtain a closed form solution for the above set or equations 

at £ — > 00. 



t(oo) 
A(oo) = (8 + 

t/(oo) „£ (8+A 2 )^(8+A 2 ) 





^ + (8+A 2 



-s/2 



s/2\ 



V(oo) = 



.n/2 



j£ 



-%/2 



_ V2[/ (8+A^)2(-(8+A^) 4 +(8+A^^-A v2 ) 



■v/2'i 



A 1 



x/2 



(18) 
(19) 

(20) 
(21) 



([4a4a> C i<?' C *+l<" C lf3 C 3 + 10} + \4a4a' C ™' C i+^ > C J+1/3 C ^] 

~i C i-lo- C io-' C «o'" :; io'' C j/3 C J + 1 /3-l — l c i-lcr C ia' c i<T' C i"" i C j+l/3 C J/3-l 

+ /l.C.). (11) 

First and third terms are irrelevant, while second and fourth 
terms respectively become 



Eq. ( fT8l shows that the Hamiltonian is diagonal in infinity. 
The hopping term has been renormalized to zero. This has be- 
come possible because of the existence of a non-zero value of 
A(£ = 0) (Eq.fPfl). Presence of any non-zero A changes the 
symmetry of the problem, and also introduces a new energy 
scales which makes the life different from that of the usual 
Mott-Hubbard physics. In case of a pure Hubbard model, one 
needs to take split the hopping terms into those which change 
the double occupancy, and those which do not 1301 bill . But 
here due to presence of a new energy scale, A, the whole hop- 
ping term can be renormalized to zero. 

Eqs. ([T9l l. (1201 . and (t2TT > indicate that the result of integrat- 
ing out the hopping term in the Hamiltonian renormalizes the 
diagonal terms. Integrating out the hopping term also induces 
new coupling which are not of the form ©, and hence irrele- 
vant. Therefore the fixed point of the one dimensional IHM is 
described by the effective Hamiltonian 

A(oo) ^-^ • r— % 

r 2^'~ 1 ) nia + U (°°) Z^ n it n U 



H(oo) 



-V(oo) 






"j(T ^Z+ 1(7/; 



(22) 



( 



4-la4 l < J >£ la ,C i -i a 5 l . j + 1 5 a p - ct CT ct r ,c i0 ./Ci <7 &_i,j<k / 3j 



+irrelevant terms. 



where ma is the number of electron in site i with spin a. Also 
note from Eq. (|2TT > that for any U, A > 0, V^(oo) is always 
negative. Hence the induced n.n. Coulomb interaction is at- 
tractive. As a consistency check we note that our effective 
Hamiltonian in the limit U, A ~^> 1 Eq. ( 1221 reduces to the 
atomic limit (t = 0) of Eq. (Q~|i. 
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FIG. 2: tonicity for A = 20, t = 1 and L = 20 sites. 
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FIG. 3: Spin and charge gaps versus U for A = 20, t — 1 and 
L = 20 sites. 



PHASE TRANSITIONS 



Spin and charge gaps 



Ionicity 

Fig. [T|represents two limiting case of the Hamiltonian ||2). 
For the band insulator, where both electrons of spin | and { 
prefer the B sites with energy —A/2, the unit cell is fully 
polarized. In the opposite limit of the Mott insulator, there 
is only one electron per site, and the unit cell polarization is 
zero. Therefore we define the ionicity for IHM as ub — ua 
where tib and tia are the density of electron in B and A sites 
respectively: 



n B = 



ua 



o-.iCB 



(23) 
(24) 



cr.iSA 



The ionicity for the IHM is depicted in Fig. [2] As can be 
seen in the figure, the ionicity is 1 for U < U Cl , indicating 
a band insulating state. For U > U C2 , the ionicity becomes 
zero, characterizing a Mott insulator. The intermediate region 
U Cl < U < U C2 is characterized by 1 > ns — Ua > 0. This 
basic physics in contained in attractive nature of the induced 
V term in Eq. (l22l . Without this term, there is only one phase 
transition at U c = A (Eq. |2). However, the renormalization 
process leading to generation of V term, opens up an inter- 
mediate state allowing for the smooth decrease of the ionicity 
from 1 (Band insulator) to (Mott insulator), by partially po- 
larizing the unit cell. 

In this picture the reduction in the degree of polarization 
of the unit cell to intermediate values facilitates the charge 
transfer between the two sites of each unit cell. Therefore the 
resulting intermediate state must be a metal. 



To investigate the intermediate state further, we numeri- 
cally calculate the charge and spin gaps for the IHM. Spin 
and charge gaps are defined as follows 



N_ 
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1 



1) 



A„ = 



E (f +!,£+ )+E (- 
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E ( 






(25) 
(26) 



where Eq{N^, N±) is the ground state energy in a sector with 
iVj {Ni) electron in spin up (down) state. Since the uni- 
tary transformations do not affect the level spacing, the spin 
and charge gaps defined above will be the same for effective 
Hamiltonian d22l as well as the original Hamiltonian (fTJi. For- 
tunately the effective Hamiltonian ( f22l is classical and hence 
free of fermionic minus sign issues. Therefore the ground 
state energy Eq(Ni,Ni) in each sector can be numerically 
calculated with straight forward algorithms. 

In Fig. Owe have plotted the spin and charge gaps versus U 
for a fixed value of A = 20. The numerical calculation is done 
for L = 20 lattice sites, subject to periodic boundary condi- 
tions. This figure shows two transition points at U Cl = 19.7 
and U C2 = 20. For U < U C1 spin and charge gap are positive 
and identical. The gaped charge excitations characterizes an 
insulating state. However, since the spin gap coincides with 
the charge gap, the resulting insulating state must be a simple 
band insulator. For U Cl < U < U C2 spin and c harge ga p are 
both zero, characteristic of metallic states H4L 1 1 0L 1 1 3L 1 1 4fl . 

Finally, for U > U C2 charge excitations will become gaped, 
while spin excitations remain gap-less. Low-energy spin exci- 
tations, with gaped charge excitation is characteristic of Mott 
insulators. This picture is in agreement with the one emerging 
from Fig. |2j according to which, rii = ^ CT ni a = 1. 

For A = 1 the transition points occur in U Cl = 0.3 and 



U C2 = 0.7 and for A 



2 in U C1 = 0.9 and U C2 = 1.7. 



The phase diagram depicting the intermediate metallic region 
is shown in Fig. |U In absence of t term, the metallic region 
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FIG. 4: Phase diagram of the IHM base on the effective Hamiltonian. 



is characterized by line U = A In this case both U Cl and 
U C2 coincide and the metallic region shrinks to a line. How- 
ever, inclusion of quantum fluctuations via t term, induces 
an attractive V term which eventually broadens the metallic 
line U = A into the narrow region shown in Fig. |4] In the 
limit U,A^>t where the effective Hamiltonian reduces to the 
atomic limit (|2), the intermediate region for a fixed A shrinks 
to zero and the metallic remains single point at U c = A 1 1 (III . 
Therefore the large A, U 3> t limit of the Fig. @] is a line on 
which band insulating and Mott insulating phases are degen- 
erate. 

In the ground state |^ ) = l^ r o( 00 )) of (I221 l. the quan- 
tum fluctuations are frozen. The ground state of the original 
fermions is given by |\&o) = U(oo)\^> ). When this transfor- 
mation is applied, it is expected to restore the low-energy spin 
spectrum of the Mott insulating state, as well as the excitation 
spectrum of the resulting metallic state. 

The question of whether the intermediate metallic state 
is a Fermi liquid, or non-Fermi liquid can be inves- 
tigated by studying the dynamic correlation functions 
(*n|Oj(£)Oo(0)|*o)- But since it is difficult to calculate |* ) 
of the original fermions, one can use the generator to set up 
flow equations for the observable O, to calculate instead the 
quantity {^ \Oj(t)Oo(0)\^o) . As far as the calculation of 
static quantities such as various energy gaps presented in this 
paper is concerned, the way one does the unitary transforma- 
tion will not be important. But for dynamic correlation func- 
tions in order to capture the incoherent features in the elec- 
tronic structure of strongly correlated systems such as IHM, 
one needs to organize the hopping processes depen ding on 
whether they change the double occupancy or not l3(Al31 l. 

To summarize, we have presented a strong coupling treat- 
ment of the ionic Hubbard model within the CUT approach. 
This enabled us to map the IHM into an effective Hamilto- 
nian containing only commuting variables rii a . Getting rid of 
fermionic minus sing problems in this way, enables us to in- 
vestigate the energetics of the model Hamiltonian, suggesting 
the phase diagram presented in Fig. [4] 
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